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O 1 Introduction 

o 

^ The study of dynamical ordering phenomena and consensus formation in ini- 

Q tially disordered populations is a central problem in the statistical physics 

' ^ approach to social and natural sciences, whether the broad idea of consensus 

is referred to opinions, voting intentions, language conventions, social habits 
or inherited genetic information. Indeed, people tend to align their opinions 
[1] , segregated populations gradually lose their genetic diversity [2] , and differ- 
ent social groups spontaneously develop their own seemingly arbitrary traits 
t-H of dress or jargon [3]. Understanding how global order can emerge in these 

situations, self-organized by purely local interactions, represents an important 
theoretical and practical problem. Our ability to grasp these issues has been 
mainly driven by the analysis of simple statistical models, which capture the 
(T) essential ingredients of a copying/invasion local dynamics. The simplest of 

those copying/invasion processes are the voter model [4] and the Moran pro- 
f*|* > ) cess [5], which focus principally on social [1] and evolutionary [6] dynamics, 

CN) respectively. 

These two models describe a population as a set of agents, each one carry- 
ing a state (opinion, trait, genome) represented by a binary variable a = ±1. 
At each time step, an ordered pair of adjacent agents (i, j) is selected at ran- 
dom. In the voter model, as a paradigm of copying processes, the system is 
updated as Oi := aj, the first agent copying the state of its neighbor. The 
voter model can thus be conceived a very simplistic model of opinion forma- 
tion in society in which individuals select their beliefs by the (admittedly not 
very realistic) procedure of just imitating one of their neighbors. On the other 
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hand, in the Moran process it is the neighboring agent the one who copies 
the state of the first agent, crj := <Tj, or, from another perspective, the state 
of the first one "invades" the neighbor in contact. The Moran process repre- 
sents thus a simple approximation to the evolutionary dynamics of a haploid 
population, constrained to have a fixed size [2]. 

In finite, initially disordered systems, and in the absence of bulk noise (i.e. 
agents spontaneously changing their state [1]), stochastic copying/invasion 
dynamics lead to a uniform state with all individuals sharing the same value 
a, the so-called consensus. Being the two final states symmetrical in prin- 
ciple, the dynamical evolution towards either state will only depend on the 
initial configuration. In order to characterize how the consensus is reached, 
two quantities are usually consider. The exit probability E(x) is defined as 
the probability that the final state corresponds to all agents in the state +1 
when starting from a homogeneous initial condition with a fraction x of agents 
in state +1. Accordingly, the consensus time T^(x) is defined as the average 
time required to reach consensus, independently of its value, in a system of 
N agents. 

Voter-like processes were originally considered in regular topologies. In this 
case, the voter model remarkably is one of the few stochastic non-equilibrium 
models that can be exactly solved in any number of dimensions [7, 8]. As it 
turns out, its symmetries imply the ensemble conservation of magnetization 
m = J2i Vi/N, which in turn implies that the exit probability takes the linear 
form E(x) = x. On the other hand, the consensus time can be seen to scale 
with the number N of agents as T N ~ iV 2 in d = 1, T N ~ Nln(N) in d = 2 
and T N ~ N for d > 2 [9]. 

The results of the voter and other copying/invasion processes in the con- 
text of social and evolutionary dynamics acquire a larger relevance when they 
are considered in systems endowed with realistic, non-trivial topologies. In 
fact, the strong heterogeneity of social and ecological substrates is better en- 
coded in terms of a complex network [10], rather than by a homogeneous 
d— dimensional lattice. Therefore, a large theoretical effort has been recently 
devoted to uncover the effects of a complex topology on the behavior of the 
voter model, as well as general dynamical processes [11, 12]. At the most basic 
level, while voter and Moran models are equivalent at the mean-field level for 
regular topologies, if the connection pattern is given by a complex network 
they behave differently, since the order in which the interacting agents are 
selected becomes relevant [13, 14]. Additionally, the heterogeneity in the con- 
nection pattern, as measured by the degree distribution P(k), defined as the 
probability that a randomly chosen agent is connected to k other agents (has 
degree k) plays a relevant role in the scaling of the consensus time, especially 
for scale-free networks [15] with a degree distribution scaling as a power-law, 
P(k) <~ fc~ 7 [13, 16, 14]. A theoretical understanding of the behavior of the 
voter model on heterogeneous networks was finally put forward by Redner and 
coworkers [17, 18, 19], who showed that for the voter model the consensus time 
scale with the system size as T N ~ N(k) 2 /(k 2 ), i.e. inversely proportional to 
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the second moment of the degree distribution (A; 2 ) = ^2 k k 2 P(k). In scale-free 
networks with degree exponent 7 < 3, the second moment (fc 2 ) diverges and 
becomes size dependent, implying thus a sublinear growth of the consensus 
time, as previously observed in numerical simulations [13, 16], and contrarily 
to the homogeneous mean-field expectation. 

When discussing the role of heterogeneity in social or evolutionary con- 
texts, a relevant question is whether the complexity of the substrate alone is 
able to encode the heterogeneity of a realistic dynamical process of social or 
environmental relevance. The common objection to extremely simple models 
is that in reality individuals behave and relate to their peers in different ways, 
i.e. they are heterogeneous both in the way in which they interact with others 
and in the way in which they react to these interaction. For example, in a 
social context, it could be the case that some agents are more reluctant to 
change their opinion (zealots), while some agents can assign different impor- 
tance to the opinion of their different neighbors (i.e. close friends can be more 
trusted than casual acquaintances). Lately, different variations of voter-like 
models have been put forward, in an effort to take into account the intrinsic 
variability of agents and their individual propensity to interact with peers 
[18, 20, 21, 22, 23]. 

In this chapter, we describe a generalization of the voter model on complex 
networks that encompasses different sources of degree-related heterogeneity 
and that is amenable to direct analytical solution by applying the standard 
methods of heterogeneous mean- field theory [11, 12]. Our formalism allows for 
a compact description of previously proposed heterogeneous voter-like models, 
and represents a basic framework within which we can rationalize the effects 
of heterogeneity in voter-like models, as well as implement novel sources of 
heterogeneity, not previously considered in the literature. 

2 A generalized heterogeneous voter-like model on 
networks 

We consider a generalized heterogeneous voter-like model, given as a stochastic 
process on networks, and defined by the following rules [24]: 

• Each vertex i has associated a given fitness /j. 

• A source vertex i is selected at random, with a probability fi/J2j fj> i- e -; 
proportional to its fitness /j. 

• A nearest neighbor j of i is then selected at random. 

• With probability Q(i,j), i copies the state of vertex j with. Otherwise, 
nothing happens. 

The fitness function fi affects the probability that the given node i is cho- 
sen to initiate the opinion-update process at a given time t [18]. In the case 
of copying dynamics, it measures the propensity of a given node to change its 
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state. Once the individual i is chosen and a neighbor j is selected, the proba- 
bility Q(i,j) measures how probable the actual update is, introducing a weight 
in the adjacency relation between the two individuals [23]. We can easily check 
that the standard voter model is recovered by setting fi = Q(i,j) — 1, while 
the Moran process corresponds to fi = 1 and Q(i,j) — ki/kj. Other variations 
of copying/invasion dynamics can be recovered by the appropriate selection 
of the fi and Q(i,j) functions. 

Recently, a generalized formalism for the class of copying/invasion voter- 
like models on networks has been proposed [25, 26] in which the process is 
identified by the copying rate dj , encoding the full structure of the contact 
network and the stochastic update rules, and that in our case takes the form 

C^ = =^^Q(i,j), (1) 

where is the adjacency matrix of the network, taking value 1 if vertices i 
and j are connected by an edge, and zero otherwise. Within this formalism, 
exact results can be obtained, but at the expense of computing the spectral 
properties of the copying rate matrix Cij , a highly non-trivial task unless the 
matrix Cij has a relatively simple form. 

Here we follow a different path, applying to our model the technique of 
heterogeneous mean-field theory, which leads to simple estimates for central 
properties such as the exit probability and the consensus time in a rather 
economical way. While this technique is known to be not exact in several 
instances, it still nevertheless able to account with a reasonable accuracy for 
the results of direct numerical simulations of the model. 

3 Heterogeneous mean-field analysis 

The analytical treatment of voter-like models on complex networks is made 
possible by the heterogeneous mean-field (HMF) approach, which has tradi- 
tionally provided a powerful analysis tool for dynamical processes on hetero- 
geneous substrates [12, 11]. Two main assumptions are made: (i) Vertices are 
grouped into degree classes, that is, all vertices in the same class share the 
same degree and the same dynamical properties; [ii) The real [quenched) net- 
work is coarse-grained into an annealed one [11], which disregards the specific 
connection pattern and postulates that the class of degree k is connected to 
the class of degree k! with conditional probability P(k'\k) [27]. In general, 
the HMF approach allows for simple analytic solutions. In the case of voter- 
like models it has proved remarkably powerful in estimating the quantities 
of interest, showing reasonable agreement with numerical simulations in real 
quenched networks [17, 18, 19, 20, 23]. 

Following the standard HMF procedure, we work with the degree-class 
average of the fitness function and microscopic copying rate. Averages are 
taken over the set of vertices with a given fixed degree, i.e. 
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~^ NP(k) ^ ^ ~~ '^ fe ' ^ 



where i € fc denotes a sum over the degree class k and P(k) is the network's 
degree distribution. Thus, represents the fitness of individuals of degree 
fc, assumed to be the same, depending only on degree, for all of them, while 
Q(k, fc') represents the probability that a vertex of degree fc copies the state 
of a vertex of degree fc'. Analogously, the contact pattern is transformed ac- 
cording to the method described in Ref . [28] , obtaining 

ay [NP(k)r 1 E l ekE 1 ek'^ ,,,,,, () 

Disregarding the microscopic details of the actual contact pattern, our gener- 
alized voter model is thus defined in terms of the mesoscopic copying rate 

C(k,k>)=-^P(k>\k)Q(k,k>). (5) 

Here and in the following me adopt the convention (•) = ^2 k P(k)(-) In 
order to provide a quantitative measure of the ordering process, we shall 
consider the time evolution of the fraction of vertices of degree fc in the state 
+1, Xk- Transition rates for x k will be given by the probability iT(fc;cr) that 
a spin in state a at a vertex of degree fc flips its value to —a [17, 19, 23] in 
a microscopic time step. It is easy to show that, from the definition of the 
generalized voter model, such probabilities can be written as 

7T(fc; +1) = x k P(k) - x k ,)C{k, fc') (6) 

fc' 

7T(fc; -1) = (1 - x k )P(k) 2 x k >C{k, fc'), (7) 

fc' 

thus leading to the rate equation [23] 

ik = Zr(fc;-l^-/7(fc;+l) ^ £ k , ){xk/ ^ (g) 

Given the very broad definition of the interaction rate Q(k, fc'), a solu- 
tion to the problem at hand cannot be easily provided in closed form, unless 
further assumptions are made. We well thus make the following additional 
assumptions: 

(i) Dynamics proceed on uncorrelated networks, i.e. [29] 

P(k'\k) = (9) 
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(ii) The interaction rate can be factorized as 

Q(k,k') = a{k)b(k')s(k,k'), (10) 

where s(k, k') is any symmetric function of k and k' . 

This simplified form encompasses a broad range of voter-like dynamical pro- 
cesses, including most previously proposed models and a variety of novel ap- 
plications, with the remarkable advantage of being promptly solvable, as it 
will become clear in the rest of this Chapter. 

In order to provide a general solution in the most compact notation pos- 
sible, we rewrite the rate equation as 



Xk 



J2P(k')r(k,k')(x k ,-x k ), (11) 



where we have defined r(k, k') = u(k)v(k')s(k, k') and 

We analyze the behavior of the linear process at hand in the canonical way, 
by first determining the inherent conservation laws [17, 18, 19]. We define a 
generic integral of motion uj[xk{t)] for Eq. (11), such that cLo/dt = 0. By 
definition of total time derivative 

dui . \ div . 

-=V xW .x = ^— x fe = 0. (13) 

k 

In analogy with previous results [17, 18, 19], we look for conserved quanti- 
ties that are linear in Xk imposing dco/dxk = Zk independent of x k , so that 
conserved quantities will be given by 

uj = z-x = ^2z k x k , (14) 

k 

where Zk is any solution of ^2 k Zkik = and x k is given by Eq. (11). It is 
easy to prove that z k oc P(k)v(k)/u(k) always satisfies the above condition, so 
that a conserved quantity is found up to multiplicative factors and additive 
constants. Upon imposing zj. = 1 as one of the possible normalization 
conditions, the conserved quantity becomes 

(v(k)/u(k) Xk) 
" = <„(*)/„(*)> • (15) 

In analogy with the simplest definition of the voter model [17] the conserved 
quantity bears all the information required to calculate the exit probability 
E, which we previously introduced as the probability that the final state 



Generalized voter-like models on heterogeneous networks 



7 



corresponds to all spins in the state +1. In the final state with all +1 spins 
we have u) = 1, while lo = is the other possible final state (all —1 spins). 
Conservation of u) implies then u = E ■ 1 + [1 — E] ■ 0, hence 

(v(k)/u(k)x k ) 

E -"- mm) • (16) 

Starting from a homogeneous initial condition, with a given density x of ran- 
domly chosen vertices in the state +1, we obtain, since w = x, 

E h (x) = x, (17) 

completely independent of the defining functions a, b, and s, and taking the 
same form as for the standard voter model [1]. On the other hand, if the 
initial condition corresponds to a single seed, that is an individual +1 spin in 
a vertex of degree fc, 

UJ - El{k) - N{v{k)/u{k)Y (18) 

which does not depend on the functional form of the symmetric interaction 
term s(k, k'). 

Eq. (11) predicts that the set of variables of x k rapidly converge to a steady 
state. It is easy to see that any choice of x k that is constant in k is a solution 
to the steady state condition x k = 0. This solution is unique and does not 
depend on initial conditions if the square matrix P(k')r(k, k') is irreducible 
and primitive (it certainly is when working with positive rates, which we will 
do in the following) [30]. If we call the steady state x°°, then it is easy to 
prove that 

u = y^^Zk'Xk' = x°°, (19) 
k> 

that is the steady state value for x k equals the conserved quantity. Such result 
is well-known in simpler formulations of the voter model and becomes crucial 
in the computation of the consensus time, even in our generalized case. 

As we noted above, the convergence to the steady state distribution occurs 
on very short time scales. As soon as the steady state is reached, stochastic 
fluctuations become relevant and the systems begins to fluctuate diffusively 
around this value, until consensus is reached in one of the two symmetric 
states. Such fluctuations are integral to finite systems and occur at long time 
scales, making this time-scale separation possible in large enough systems. 
In the light of such considerations, the average consensus time Tjv(x) for a 
system in a generic steady state x can be derived extending the well known 
recursive method to our general case [19]. At a given time t, Tjv(x) must equal 
the average consensus time at time t + At plus the elapsed time At = 1 /N 
that is, in our notation, 

Tjv(x) = TTTjv(x) + n(k; s)Tjv(x + Ax'*') + At, (20) 

k,s 
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where 5=1 — J2k s H{k; s) is the probability that no state change occurs, 
while the sum is the weighted average over possible state-updates x — > x + 
Ax' 1 '. The variation Ax' 1 ' is a vector whose all components are zero except 
for the k-th, which equals the update-unit A k = [NP(k)]^ 1 . Expanding to 
second order in A k , taking x k — lo as the initial state and changing variables 
such that d/dxk — Zkd/duj we obtain the backward Kolmogorov equation 

z T Tz . ,d 2 T N 

leading to 

Tat = -N cS [u\nuj+ (1 - w) ln(l - u)] (22) 

where we have defined the effective system size iV e ff = ^/Sfe fe' z fe^(fc, k')zk>, 
which, in the case of generalized voter dynamics, Eq. (12), becomes 



(/(*)>(*> (/S& j) 
's{k,k')kb{k)^ 

where we have defined ((•)) = £ fcfe , P(k)P(k'){-). 



2 



N eB = N / , , (23) 



4 Particular cases 

With the formalism developed above, we can easily recover several of the 
variations of the voter model proposed in the past. Let us look at some of 
them in following section. 

4.1 Standard voter model and Moran process 

The standard voter model and Moran process can be recovered by setting 
a{k) = b(k) = f(k) = s(k,k') = 1 and a(k) = k, b(k) = fc" 1 and f(k) = 
s(k,k') = 1, respectively. In this case the known results [19, 17] are recovered. 
Thus, for the voter model, the conserved quantity is 

- = £*^P^), (24) 

the exit probability starting from a single +1 vertex is 

Mk) = -±r y (25) 

and the consensus time takes the form 

(k) 2 

T N (u) = -Nj^r [w Inw + (1 - w) ln(l - w)] . (26) 
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On the other hand, for the Moran process we have 



(fc-i) 



1 




(27) 



Ei(fc) 



1 1 



(28) 



T n (uj) 



kNik- 1 ) 

-N(k) (k-^x [u In(w) + (1 - u) ln(l - w)] . 



(29) 



The difference between voter and Moran dynamics is quite evident here. By 
looking the expressions for the conserved quantities, Eq. (24) states that for 
the voter model densities are weighted with a factor k/(k) which compensates 
the tendency of small degree nodes to change their state, whereas in the Moran 
process the exact opposite occurs, being the density in Eq. (27) balanced by 
the factor fc _1 /(fc _1 ). This fact translates in the different forms of the exit 
probability, Eqs. (25) and (28). For the voter model, a single mutant opinion 
can spread to the whole system more easily when it first starts in a vertex of 
large degree, due to the fact that large degree vertices are copied from with 
larger probability. On the other hand, a mutant in the Moran process is able 
to spread faster if it starts on a low degree vertex, owing to the corresponding 
fact that low degree vertices are invaded with low probability. 

4.2 Voter model on weighted networks 

The extension of voter model to weighted networks [23] is motivated by those 
situations in which the strength of a relation can play a role in the process of 
opinion formation. In this sense, weights would reflect the evidence that the 
opinion of a given individual can be more easily influenced by a close friend 
rather than by a casual acquaintance. In a weighted network, the voter model 
is defined as follows: At each time step a vertex i is selected randomly with 
uniform probability; then one among the nearest neighbors of i, namely j, 
is chosen with a probability proportional to the weight u>ij > of the edge 
joining i and j. That is, the probability of choosing the neighbor j is 



Vertex i is finally updated by copying the state of vertex j. If the weights de- 
pend on the degree of the edge's endpoints, Wij — g(ki,kj), with g(k,k') 
a symmetric multiplicative function, i.e. g(k,k') = g s (k)g s (k') [23], voter 
dynamics is recovered by setting s(k, k') — a(k) — f(k) — 1 and b(k) = 
g s {k)(k) j (kg s (k)) , which leads to an invasion exit probability 



E\(k) = 



kg s (k) 



(31) 



N(kg s (k)) 



and a consensus time 
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TjvM = ~ N W^W) [uj1ulj + {1 - w) ln(1 - w)] ' (32) 



with the conserved quantity 

k'g s (k')P(k') 



,, (kg a (k)) 



x k ,(t). (33) 



In order to provide an example of weighted voter dynamics, we can consider 
the special case of weights scaling as a power law of the degree, g s {k) — k 8 
on a scale-free network with degree distribution of the form P{k) ~ fc~ 7 . The 
consensus time starting from homogeneous initial conditions, Xk(0) = 1/2 
takes the form 

T N (l/2)=Nln(2))^- y (34) 

From this expression, we can obtain different scalings with the network size N, 
depending on the characteristic exponents 7 and 9. Considering only 7 > 2 
and using the scaling behavior of the network upper cutoff k c <~ N 1 / 2 for 
7 < 3 and k c ~ iV 1 /^" 1 ) for 7 > 3 [31], we are led to different regions of 
behavior for TV (1/2): 

1. If 9 > 7 - 2 both (k 1+e ) and (fc 2 + 20 ) diverge. In particular, (fc 1+e ) ~ 
k 2 + B -« and (k 2 + 29 ) ~ fc c 3 + 2 ^. Thus 

Tat ~ iVA^. (35) 

If 7 < 3, k c ~ TV 1 / 2 , and TV - JVO-t)/ 2 . If 7 > 3, then fc c - ArV(7-i) ; 
and Tjv ~ const. 

2. If 7-2 > 9 > (7-3)/2, then (/c 1+e ) converges and (fc 2 + 29 ) diverges. Thus 

TV ~ Nk2~ 20 - 3 . (36) 
If 7 < 3, Tjv ~ iV(7-2e-i)/2. if 7 > 3, then Tjv - 7V 2 (7-e-2)/( 7 -i)_ 

3. If 6< < (7 - 3)/2, then both (fc 1+e ) and (fc 2+2e ) converge, and we have 

TV ~ N. (37) 

These scaling relations contain several interesting aspects. Among these, 
it is worth highlighting that, in region 1, and for 7 > 3, the analytical results 
predict a constant scaling for the consensus time, which therefore docs not 
depend on the population size. As a consequence, in the thermodynamic limit 
the ordering process is instantaneous. However, this turns out to be true only 
at the mean field level, i.e. on annealed networks. Simulating the process on 
quenched networks produces in fact different results, with a clear dependence 
of the consensus time on N [23]. Henceforth, this is a typical case showing 
the limits of the mean field approach in predicting the behavior of phenomena 
occurring on true finite networks with fixed connections and disorder. Another 
interesting feature concerns the special value 7 = 3. When 9 > 0, this value 
separates distinct scaling behavior, while as 9 < it ceases to be a frontier, 
the scaling of the consensus being linear in TV on both sides of it. 
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5 A practical example: variable opinion strengths 

The strength of the proposed generalized formalism resides in the possibility 
of deriving simple HMF solutions to new variations of the standard voter-like 
models, which were not studied in the past and whose solution would have 
been too hard to compute with more exact techniques. As an application of 
our formalism, we consider the case of voter dynamics in a society in which 
certain individuals are more likely to align their opinions with those of their 
neighbors than others [18]. In particular, the propensity of a certain individual 
to change opinion will depend on the strength of his/her social ties, that is 
on his/her number of neighbors. We can accomplish such a description in 
our formalism, by encoding this dependence in the fitness function f(k) and 
assuming a{k) = b(k) = s(k, k') = 1 for the sake of simplicity. The HMF 
solution to such a problem will then be easy to derive. Following the steps 
illustrated in the previous sections, the conserved quantity reads 

the exit probability for an individual seed in state +1 

_ k/f(k) 

E ' nWJWY (39) 

and the consensus time starting from a homogeneous state uj 

Tn(co) = -N { -l^W^ [u lnu, + (1 - W ) ln(l - uj)]. (40) 

We can provide a practical example by choosing f(k) — k a . If a < 0, 
less connected individuals are more likely to change opinion, and connected- 
ness can be interpreted as a measure of social self-assurance. If a > 0, more 
connected individuals appear more vulnerable to opinion variability. Inter- 
estingly, the standard voter model is recovered for a = 0, where no agent 
heterogeneity is postulated. The conserved quantity assumes the simple form 
uj = (k 1 ~ a Xk)/(k 1 ~ a ) and the consensus time starting from an initial condi- 
tion uj = 1/2 is then 

/h.a\ (h.\-a\2 

T N {l/2) = NH2) { ^ 2 _ a) } . (41) 

The size scaling behavior of the consensus time can be derived from Eq. (41), 
following the procedure devised in the previous section for the voter model 
on weighted networks. The results are illustrated in compact form in Fig. 1, 
where a phase diagram for the variables 7 and a is shown. The phase diagram 
conveys great insight into the role played by the parameter a. Focusing on 
the region with 7 > 2, which is of greater interest in the study of dynamical 
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processes on complex networks, we obseve that for 3—7 < a < 7— 1 the simple 
scaling relation Tjv ~ N is recovered. In analogy with a similar phenomenon 
observed in the case of weighted networks in Section 4.2, the value 7 = 3 
does not act as frontier and the same scaling law is observed, regardless of 
the value of 7. Larger values of a break the balance that ensues the Tjy ~ N 
behavior and lead to non-trivial size dependence of the consensus time, with 
exponents that depend on the degree distribution exponent 7. If we restrict 
our analysis to the case of scale-free networks, corresponding to the 2 < 7 < 3 
strip in the phase diagram, we can easly see from the results in Figure 1 that 
values of a in the range \a\ < 7 — 1 lead to either linear or sub-linear size 
scaling of Tjv, whereas for \a\ > 7— 1 super-linear scaling is encountered. This 
translates into the simple observation that in scale-free networks, large degree 
selectivity makes consensus harder to reach in larger systems, regardless of the 
sign of a, i.e., of whether the degree selectivity makes low-degree individuals 
or high-degree individuals more vulnerable to opinion change. 

In order to corroborate our predictions, we have performed numerical sim- 
ulations of the voter dynamics at hand in complex networks with power-law 
distributed degrees, generated with the Uncorrelated Model (UCM) [32]. We 
check the scaling behavior of the consensus time, Tjv ~ N 13 , for three pairs of 
(7, a) values of the phase diagram, located in the regions where super-linear, 
linear and sub-linear scaling of Tjy are observed respectively. 



Generalized voter-like models on heterogeneous networks 



13 




Fig. 2. Scaling behavior ol the consensus time Tjv as a lunction of the network size 
N for different values of 7 and a. 

In Fig. 2 we compare the scaling exponent predicted by Eq. (41), /?hmf, 
with the exponent (3 q obtained by fitting the numerical simulations run on 
quenched networks. The results are summarized in the same Fig. 2. We note 
that the HMF theory predicts qualitatively well the behavior of the voter 
dynamics, in the sense that the scaling of Tjv with N is super-linear for 7 = 2 
and a = 2, linear for 7 = 2.5 and a — 1, and sub- linear for 7 = 3 and a = — 1, 
as expected. As for the exact values of the j3 exponents, the HMF prediction 
is quantitatively accurate only is certain regions of the phase diagram, where 
the annealed-network approximation appears to hold. Elsewhere, quenched- 
network effects take over and sensible deviations with respect to the HMF 
value of ft are encountered in simulations, in analogy with results for the 
weighted voter model, as discussed in Section 4.2. 

6 Conclusions 

In this Chapter, we have presented a generalized model of consensus forma- 
tion, which is able to encompass all previous formulations of copy/invasion 
processes inspired by variations on the voter model and the Moran process. 
We considered the implementation of such generalized dynamics on a het- 
erogeneous contact pattern, represented by a complex network, and derived 
the theoretical predictions for the relevant dynamical quantities, within the 
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assumptions of the heterogeneous mean-field theory. We provided a brief re- 
view of previous results that can be recovered by our generalized formalism, 
and finally we considered a novel application to the case of opinion forma- 
tion in a social network. In particular, we addressed the case in which the 
opinion strength of an individual is related to his/her degree centrality in the 
network. We found that in scale-free networks strong selectivity rules (which 
make less connected individuals much proner to change their opinions than 
more-connected ones or vice versa) lead to a steeper growth of consensus time 
with the system size, making the ordering process slower in general. Numer- 
ical simulations on quenched networks show that the HMF theory is able to 
predict such behavior with reasonable accuracy. Slight deviations from the 
theoretical predictions are encountered in certain regions of the phase dia- 
gram, but they are due to quenched-network effects that the HMF theory is 
not be able to capture. 
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